AD-AO02  194  NAVAL  RESEARCH  LAB  WASHINGTON  DC  F/G  20/9 

RADIATION  FROM  A  LOCALIZED  LANGMUIR  OSCILLATION  IN  A  UNIFORMLY  — ETC(U) 
JAN  80  H  P  FREUND*  K  PAPADOPOULOS  NASA-W-14365 

UNCLASSIFIED  NRL-HR-4153 _ SB  IE  -AO-EOOO  377 _ NL _ 


rfP*£ooo  377  t 

NRL  Memorandum  Report  4153 


Radiation  from  a  Localized  Langmuir  Oscillation 
in  a  Uniformly  Magnetized  Plasma 

H.  P.  Freund 

Science  Applications,  Inc., 

McLean,  Virginia  22101 

AND 

K.  Papadopoulos 
Plasma  Physics  Division 


January  29,  1980 


DTIC 

SELECTEI 

MAR  2  5  1980 1  | 

B 


NAVAL  RESEARCH  LABORATORY 
Washington,  D.C. 


Approved  for  public  release;  distribution  unlimited. 


SECURITY  CLASSIFICATION  of  this  RAGE  (When  Dim  gn tatad) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

S.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (mod  Sobllllo) 

RADIATION  FROM  A  LOCALIZED  LANGMUIR 
OSCILLATION  IN  A  UNIFORMLY  MAGNETIZED  PLASMA 

S.  TYRE  OF  REPORT  4  PERIOD  COVERED 

Interim  report  on  a  continuing 

NRL  problem. 

••  PERFORMING  ORO.  REPORT  NUMBER 

7.  AUTHORS 

H.  P.  Freund*  and  K.  Papadopoulosf 

•  CONTRACT  OR  GRANT  NUMBERS 

s.  perporming  organization  name  ano  aooness 

Naval  Research  Laboratory 

Washington,  DC  20375 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

NRL  Problem  67-0884-0-0 

NASA  Contract  W  14365 

II.  CONTNOLLING  OFFICE  NAME  ANO  AOONESS 

Office  of  Naval  Research,  Arlington,  VA  22217  and 

National  Aeronautics  and  Space  Administration,  Washington, 

DC  20546 

gMggrnjjrgi . .  ■  .  .  -  .  _ 

19.  NUMBER  OF  PAGES 

30 

14.  MONITORING  AGENCY  NAME  4  AOORESSfTl  dl  Ho  tool  Iron  Controlling  0111  cm) 

is.  SECURITY  CLASS,  (ol  thla  r apotl) 

UNCLASSIFIED 

nf.mB.iuijr-rwmi  111111111  m  1 

16.  DISTRIBUTION  STATEMENT  fol  Iblo  Rmyort) 

Approved  for  public  release;  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ot  fh*  abstract  antarad  In  Block  20,  II  dltlarant  from  Raport) 

IS.  SUPPLEMENTARY  NOTES 

♦Science  Applications,  Inc.,  McLean,  Virginia  22101 
fUniversity  of  Maryland,  College  Park,  Maryland  20742 

(Continues) 

It.  KEY  WORDS  (Continue  on  ravaraa  old*  II  n«e»M«y  and  Idantlty  by  block  nuotar) 

Soliton  radiation  Plasma  frequency  radiation 

Magnetized  plasma 

Strongly  turbulent  radiation  mechanisms 
Langmuir  solitons 

20.  ABSTRACT  (Conllnvo  on  nwn  old,  II  nciiuv  md  Identity  by  block  mmbot) 

/---The  radiation  at  the  first  and  second  harmonics  of  the  electron  plasma  frequency  from  a 
localized  Langmuir  perturbation  is  computed  for  the  case  of  a  uniformly  magnetized  plasma.  It 
Is  assumed  that  the  localized  perturbations  in  both  the  electrostatic  field  and  plasma  density  have 
cylindrical  symmetry  about  the  direction  of  the  ambient  magnetic  field.  The  analysis  is  initially 
performed  for  such  an  arbitrary  localized  perturbation,  and  then  applied  to  treat  the  case  of  a 
one-dimensional  Langmuir  soliton  propagating  in  the  direction  of  the  magnetic  field.  An  extensive 

(Continues)  -2 


DD  1473  COITION  Of  I  NOV  «»  IS  OBSOLETE 

S/N  010  2- 014*  6601  i  - - -  -  — — - — 

6ICUNITY  CLASSIFICATION  06  THIS  0*01  (BSm  Bah  tnlorod) 


SCCUKITV  CLASSIFICATION  OF  THIS  F»Ct  (Whtr,  D. la  En|.,.< 


18.  Supplementary  Notes  (Continued) 


HiU  research  was  sponsored  partially  by  the  Office  of  Naval  Research,  and  by  the  National 
Aeronautics  and  Space  Administration  under  contract  W-14365. 

20.  Abstract  (Continued) 


numerical  study  of  the  angular  dependence  of  the  radiation  spectrum  as  a  function  of  the  ratio  of 
the  electron  plasma  and  cyclotron  frequencies  is  described. 


SSCUHlTy  CLASSIFICATION  of  this  PtkOt(Wh*n  Of  tut •«<) 


CONTENTS 


I.  INTRODUCTION  .  1 

II.  THE  EMISSIVITY  .  2 

III.  THE  CASE  OF  ONE-DIMENSIONAL  SOLITONS  .  9 

IV.  SUMMARY  AND  DISCUSSION .  15 

ACKNOWLEDGMENTS .  17 

APPENDIX:  THE  DISPERSION  TENSOR  AND  SOURCE 

CURRENT .  18 

REFERENCES  .  21 


ACCESSION  (Of _ 

NT  18  Whit*  StcttM 

DOC  Butt  SKtiM 

UNANNOUNCED 

JUSTIFICATION  _ 


BIT 

BBWlON/WinUKin  CODES 

lit:  AyfflL  ind/or  SPECIAL 
- 1 — - 


□  □ 


I.  INTRODUCTION 


The  question  of  strong  Langmuir  turbulence  in  magnetized  plasmas  is 

important  in  studies  of  beam-plasma  interactions  in  both  space  and 

laboratory  plasmas.  Heretofore,  studies  of  strong  turbulence  theory  in 

magnetized  plasmas  have  centered  on  the  dynamics  of  collapse  and  the 

shape  and  stability  of  the  localized  structures  which  result.  However, 

such  structures  are  expected  to  have  electromagnetic  signatures  at 

harmonics  of  the  electron  plasma  frequency.  While  the  electromagnetic 

radiation  from  Langmuir  solitons  has  been  extensively  studied  for 

field-free  plasmas,1-4  there  has  been  scant  treatment  of  the  problem 

of  the  strongly  turbulent  radiation  process  in  magnetized  plasmas. 

This  problem  is  of  particular  relevance  due  to  increased  interest  in 

experimental  studies  of  electron  beam  driven  strong  turbulence  in  the 
5-9 

laboratory.  It  is  our  intention  to  address  this  question  in  the 
present  work,  and  to  derive  expressions  for  the  radiation  emissivity 
from  spiky  Langmuir  turbulence  at  the  first  and  second  harmonics  of 
the  electron  plasma  frequency. 

The  organization  of  the  paper  is  as  follows.  In  Sec.  II,  we 

derive  an  expression  for  the  emissivity  from  an  arbitrary,  cylindri- 

cally  symmetric  soliton  at  frequencies  w  =  tog,  2we  (where  a>e  denotes 

the  electron  plasma  frequency).  It  should  be  noted  that  the  treatment 

of  emission  at  the  electron  plasma  frequency  is  restricted  to  the 

limit  in  which  the  radiation  wavelength  is  much  less  than  the  scale 

length  of  the  soliton.  In  order  to  investigate  simplified  scaling 

laws  between  the  radiation  emissivity  and  the  soliton  amplitude,  we 

consider  the  limiting  case  of  one-dimensional  Langmuir  solitons  in 
Note:  Manuscript  submitted  December  7,  1979. 
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Sec.  III.  A  numerical  study  of  the  angular  dependence  of  the 
emissivity  is  also  presented  in  this  section.  In  particular,  we 


investigate  the  variation  of  the  radiation  pattern  with  wg/Qe  (where 

n  is  the  electron  cyclotron  frequency) .  A  summary  and  discussion 
e 

appears  in  Sec.  IV,  and  the  derivation  of  the  plasma  dispersion  tensor 
and  the  r  elation  source  current  is  given  in  an  appendix. 

II.  THE  EMISSIVITY 

We  assume  a  localized  Langmuir  perturbation  of  the  form 

jg(£,t)  =  V$(r,z)sina>et,  (1) 

where  V«(i(v,z)  defines  the  soliton  envelope,  and  the  ambient  magnetic 

field  B  (=B  e  )  defines  the  z-axis.  The  interaction  between  the 
0,0  o  o,z 

electrostatic  field  and  the  associated  slow  time  scale  oscillation  in 
the  plasma  density  is  implicitly  included,  and  we  use  Sn(r,z)  to  denote 
the  density  caviton.  Both  <p(r,z)  and  5n(r,z)  possess  cylindrical  sym¬ 
metry  about  the  z-axis,  and  it  is  assumed  that  i|>(r,z)  and  6n(r,z)  are 
odd  and  even  functions  of  z  respectively. 

The  radiated  power  is  defined  to  be 

T/2 

P  =  -  lira  ijr  f  dt  /d3  x  6E(x,t)  *  6J  (x,t),  (2) 

T-xo  T  J  J  *  %  % 

-T/2 

where  6E(x,t)  is  the  radiation  electric  field,  and  6J  (x,t)  is  the 

•\j  %  ^ 

source  current  due  to  the  localized  Langmuir  perturbation.  Equation 
(2)  can  be  expressed  in  terms  of  the  Fourier  amplitudes  of  6E(x,t)  and 
6J  (x,t)  in  the  following  manner 

f\jS  *\f 
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* 


■  VI"" 


00 


(3) 


P  -  -  (2ir)4  lim  ^  /d3k  f  du  6E(k,u>)  •  6J  *(k,to), 
T-w>  1  J  'Lx,  ^  ^  %s  ^ 


where  the  asterisk  (*)  denotes  the  complex  conjugate,  and  the  Fourier 
transform  is  defined  as  follows 


f(k,w)  -  (2n) 
%  % 


-4 


/d’  x/ 


dt  exp(iiot-ik»x)f  (x,t) . 
^ 


(4) 


A  self-consistent  relation  between  the  radiation  field  and  source 
current  is  derived  in  the  Appendix,  and  is  of  the  form 


A(k,a>)  •  SE(k,io) 

K  <1,  0.0/ 


—  5J  (k,w) , 

10  O/S  0/ 


(5) 


where  the  dispersion  tensor  is  given  by  (k  =  k  e  +k  e) 
r  °  o  x  ox  z  O/Z 


A(k,w)  *  (k  k  -  k2  I)  +  e(k,u>). 

fa  0/  (i)^-  %  fa  W  <v 


(6) 


In  Eq.  (6),  ^  is  the  unit  dyadic,  ^(k.to)  is  the  plasma  dielectric  ten¬ 
sor,  and  we  have  that  -  e„„  =  £,,  =  -e„„  =  ie2,  ezz  =  e3. 


xx  yy  l  xy  yx 

exz  *  ezx  *  eyz  *  Gzy  "  °*  where  el  2  1  ~  <?/(“)2  “  fie2)’  e2  2  “e  V 


to (u)2  -  ^e2)»  and  e  *  1  -  to'2  /to2.  Here  to'  =  ioe(l  +  ^nex/nD)  »  where 

n  is  the  ambient  electron  density  and  fin  denotes  the  extremum  of  fin. 
o  'ex 

This  dispersion  tensor  includes  the  strongly  turbulent  modifications  to 
the  cold  plasma  approximation,  and  has  been  derived  under  the 
restriction  that  the  wavelength  of  the  radiation  be  much  less  than  the 
scale  length  of  the  local  perturbation. 


Inverting  Eq.  (5)  to  find  <5E(k,cu)  as  a  function  of  6J  (k,io),  we 

i\,  'X,  'X.S  %  ' 


obtain 


10 


3 


(7) 


<5E(k,a>) 


.  .  A  (k,ai) 
4iri  ss\’ 


2  (£>w)  '  6^s(^“>> 


where  A(k,u>)  is  the  determinant  of  A(k,m),  A  (k,w)  is  the  trace  of  the 
%  ~  ss 

classical  adjoint  of  A(k,m),  and  a(k,w)  =  6E(k,a>)/ I 6E(k,w) I  is  the  unit 

^  %  'X,  'V.'V, 

polarization  vector.  The  power  radiated  per  unit  solid  angle  subtended 
by  k  is  found  by  substitution  of  (7)  into  (3),  and  it  can  be  shown  that 


§-  •  <2.)s  u„  i  e  r  *2  y*’u) 

d\  ^  Jo“ 


k2=k.  2 


x  *  5JS  (k,w)  |  6(k2  -  k^2),  (8) 


where ^  denotes  a  sum  over  the  wave  modes  of  the  system, 

± 

dn.  -  2 ir  sin  0d0  (where  0  denotes  the  angle  between  k  and  B  ),  and  we 

K.  r\j  %0 

% 

have  summed  over  the  contributions  of  positive  and  negative  w.  The 
appropriate  wave  modes  are  described  by 


c2k_j_2 


2ot2 (1  -  a2) 

2(l-a2)  -  S2(sin20?p)  * 


where  the  "  +  "  and 


denote  the  ordinary  and  extraordinary  modes 


respectively,  a2  =  u)'2/w2,  02  =  fieZ/w2,  and  p2  =  sin40  +  4(w/fte)2 
X(1  -  ct2)2cos20.  This  corresponds  to  the  well-known  Appleton-Hartree 
dispersion  relation  in  which  cu'  has  been  substituted  for  (oe  to  des¬ 
cribe  the  strong  turbulence  effect. 

In  this  paper  we  treat  emission  at  a>  =  u  ,  2ue>  and  write  the 
source  current  as  the  sum  of  contributions  (see  Appendix) 


4 


6J  (k,oi)  =  fij;  (k)  6(o>-u  )  +  6J'  (k)  6(w-  2u>  )  where 

>\jS  ^  e  e 


6J^(k)  =  -  (2ir)  4  -A  /d3x  exp(-ik*x)  ~ 

o 


( 


a  2 

e 


u>  ft 

e  e 


V*  +  iT^  V  +  1  Sz*7*)’  (10> 

e  e  e  e 


6jf2)  (k)  =  i(2TT)"4  ~7T  Jd\  exp(-ik-x) 


<\,S  <v 


8m  u 
e  e 


e  e 


a  ^ 


*  H  '2* +  V*J  l’* +  V 


e  e 


a)  ft  ^ 

+  i  -2-o' "2  [3I-2g  (k,2uj  )  ]  •  (e„  x  V<J>) 
a)  *-ft  ^  L  ^  V  e  J  %z 
e  e 


> 


f + AvX( 


+  ^  %'2“e}3  T  + 


fi 2 

+  - 7--- -7  (7 

0)  Z-Q  z  v 
e  e 


,»■)]. 


(11) 


and  V.  =  V  -  e  (3/3z).  If  we  write  the  square  of  the  delta  function  as 
J-  <\,Z 

62(o>-w  )  =  lim(T/2"r)6(a)-u)  )  and  evaluate  3A/3k2  for  the  appropriate  mode, 
0  T-*«  0 

then  the  expression  for  the  radiated  power  becomes 

u>2X  (k,u>) 


41 

da 


-  (2»)5E /  4»  "cZn*pc'  S(k2'k2±) 


e  2 


/ 1  \ *  ,2 

*  »|a,.,(k,u)  •  6JV  '  (k)  6(u)-u>  )  + 

^(i)  X  <V8  %  e 


[' 


5 


ft  vn* 


+  |at(k,W)  •  SJs‘ 


(k)  |  6  (oj-2o) 


In  the  evaluation  of  the  radiated  power  at  w  -  u  and  2w  which 

e  e 

follows,  we  make  use  of  the  expression10  a  =  -  i(X  X  )""2(X  ,X  ,X  ) 

^  ss  22  12  22  32 

to  describe  the  unit  polarization  vector,  where  X  is  the 

classical  adjoint  of  A.  It  should  be  noted,  however,  that  this  approach 

is  invalid  when  the  eigenvalues  of  A  are  degenerate  (i.e.,  when 

RJ 

k  2  =  k  2).  Since  this  occurs  in  the  limit  in  which  SI  -*■  0,  care  must 
+  -  e 

be  exercised  in  order  to  treat  the  field-free  case. 

A.  Emission  at  u>  -  w  : 

e 

In  this  frequency  regime  we  are  restricted  to  consideration  of 
waves  whose  wavelength  is  much  less  than  the  scale  length  of  the  per¬ 
turbation,  which  is  equivalent,  in  practice,  to  the  condition  that 

c2k,2/oo  2  >|fin  /n  |.  The  emissivity  n  (s  w  ^PCuO/dft  )  is  defined  to 
±  e  ex  o  k 

'Xj 

be  the  power  radiated  per  unit  frequency  per  unit  solid  angle  subtended 


by  k.  In  computing  rj(o)  ,8),  we  retain  only  the  contribution  due  to  the 

oscillatory  current  at  w  (i.e.,  6J^ ^  (k))  and  find 
J  e  '  %s  ~ 


n(V9)  =  64^ 


we3  /Sn 


N±(p±sin29) 

pcos^e 


N±2si 


sin0cos0  Ii 


oj  2  /  6n  J2  2  ' 

u, . 2  (N±2sin20  +  (sin29+p) 

e  e  \  o  e  > 


where 


co  no 

\  =  J"  izj'  drr  cosCk^z  cos©)  J^k^rsinO)  ^  Vj|$, 


-CD  o 


6 


(15) 


CO  CO 

tA  *  f  dz  f  drr  sin(k±z  cos9)J  (k±rsin0)  ■—  V^4>  , 
J — °°  J  o  1  o 


p2  =  sin40  +  4(we2/fte2)  (5nex/no)2cos28,  N(=ck/ug)  is  the  index  of 


refraction,  and 


n±  =  1  ~  6^ 


6n  /  6n  \ 

— —  (  1  +  — —  J 

no  \  no  / 


fT7 

ex  ,  e 


(16) 


r*  +  2^r  (sin^e^p) 


B.  Emission  at  to  -  2to  : 

e 

In  this  frequency  regime,  the  strongly  turbulent  contributions  to 
the  dielectric  properties  of  the  plasma  can  be  ignored,  and  the  Appleton- 
Hartree  dispersion  relation11  can  be  employed.  After  retaining  only  the 
oscillatory  source  current  at  2<oe  in  (12),  we  find  that 
g2  N±(p±sin20) 

n(2toe,0)  36Trme2c;Jioe  pcosz0  ^A||.li  I]|,|| 


+  A||,x  +  Ax,[[  Vll  +  Ax,x  ^ 


(17) 


where 


CO  CO 


Ijj  ij  *  J~  dz  J  drr  cos(k±zcos0)  JQ(k±rsin0)  (V$>)||2  , 

00  00 

I  =  f  dz  f  dr  sin(k±zcos8)  Jo(k±rsin0)  (V<|>) jj  (V4>)J 

IL-1  1  11 


1^  i  =  J  dz J  drrsiruk^zcosQ)  J^(k±rsin9)  (V$)j'(V$) 


— 00  o 


CO  OO 

IA  A  =  J'  dz  J~  drrcos (k^zcosS)  Jo(k±rsin8)  (V<f>)^2 


y  k^sinS^ 


rQ 

2%2cos20  ”  (U  2-q  2)(4U  2-n  ?y(Ni2sln26 


e  e  e  e 


,  (7oi  2-n  2)n  2 

.  1 _ e _ e _ e _  2fl. 

4  (u,  2-fi  2)  (4oo  2-0  2)  CSln  B+ 
e  e  e  e 


in2  0^p  , 


All  =  — -  t  N,  2sin8cos0 
l,x  w  z  ± 
e  e 


\,ii =  u  k±cos9  (N±2sin2e  +  (N±2sin2e  - 

e  e  \  e  e 


n  2  6oj  2-n  2  > 

e  e  e  t 


/.  610^ 

^x,x  "  2  k±sine(2  N±2cos20  -  To)  2-i"2)'(4w''^n_2y(N±2si 

\  e  e  e  e 


(5u  2-n  2)fi  2  \ 

+  4  (o>  2)  (4o)  2-0  2)  (Sin2e?P})  * 

e  e  e  e  / 


In  addition,  p2  =  sin^S  +  9(u  2/fi  2)cos20,  and 

e  e 


2  i  - 

n±  "  1  12w  2  -  2 n  2(sin28^p) 


III.  THE  CASE  OF  ONE-DIMENSIONAL  SOLITONS 


We  choose  to  apply  the  expressions  derived  in  Sec.  II  to  the  case 
of  one-dimensional  solitons  since  this  is  the  only  regime  which  is 
analytically  accessible,  and  for  which  the  angular  spectrum  of  the 
radiation  may  be  obtained  in  a  relatively  straightforward  manner.  We 
assume  the  electric  field  and  plasma  density  perturbations  to  be  of  the 
form12 


V<j>(r,z)  =  E(r)  sechLic(r)z]e  , 

'vz 


(21) 


and 


(r .? g.)  =  _  g  K2(r)A  2  sech2Lic(r)z], 


(22) 


where  A^  is  the  electron  Debye  length,  l/ic(r)  characterizes  the  parallel 
scale  length  of  the  perturbation, 


o  e  \  e  ' 


(23) 


W(r)  —  E2(r)/8ir,  and  are  the  ion  temperature  and  ratio  of  specific 
heats,  and  Te  is  the  electron  temperature.  In  addition,  we  assume  that 
<(r)  =  KQexp(-  r/p),  where  P (»  kq  J)  characterizes  the  transverse 
scale  length  of  the  perturbation.  It  must  be  remarked  that  in  order  to 


neglect  the  radial  component  of  the  soliton  field,  we  must  have  that 

2 


II  o  a 


(24) 


and 


Jo(kKo_1)l7/l  ?>  tan(k<o_1)J  (kKo_1)  ^7  ?  lYl* 

ii  1  e  e 


(25) 


9 


...»  Alli- 


r 


As  a  result,  we  cannot  treat  the  case  in  which  is  arbitrarily  close 

to  ft  by  means  of  Eqs.  (21)  and  (22). 
e 

Use  of  (21)-(23)  immediately  yields  the  following  expressions  for 
the  emissivity  at  and  2u)^ 


_  2  /v  \3  n  T  /6n  \  -2  /  y.T.  \-1 

■«V>-£-(r)^T  (nf  (-t1) 


sin2  9 


^  N±5(p±sin20) 

*  P  Tl±2 


(26) 


and 


n(2ue,e) 


2  /v  \3  n  T  l 

=  f-  (/-)  -p  tan2ej} 


kN ,  (p±sin20) 

=— r - V  2(e)I  2, 

D  ±  2± 


(27) 


where  v  2  =  T  /m  , 
e  e  e 


3  2w  4  - 

V9)  *  I  N±2“s29  -  (V/-ae^)(4,J^-a^)<11±2sln29  -  4> 


7u>  2  ft  2 

+  4^e^)(4me^-ne^)  (sin20^>’ 


(28) 


and  the  source  integrals  are 

2  2 


/■  /  k  2cos20\  .  .  / irk,  cos  0  \ 

drr  ^1  +  j  rr  J0(k±rsine)sech^  ^ (r) "  J  , 

/k  2cos0  .  /irk^cos©^ 

drr  'V(7T  ;rr  vvsin0)csch  • 


(29) 


(30) 
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The  only  regime  in  which  (29)  and  (30)  can  be  integrated 


analytically  is  the  case  for  which  k±  <  kq.  However,  it  is  important  to 

recognize  that  n(u>e,0)  has  been  derived  subject  to  the  condition  that 

k,  >  k  (which  is  equivalent  to  6  v2  N  2  >  c2  1 6n  /n  I )  and  no  further 
x  o  e  ^  ex  o 

analytic  reduction  is  possible  for  the  case  of  emission  at  u>e  (29) . 
Turning,  therefore,  to  the  case  of  emission  at  2wg,  we  find  that  when 


6 


<< 


(31) 


the  emissivity  becomes 


n(2ve)  -  ^1^ 


✓v  \4  c3W  /  y.T.V 


(32) 


where  Wq  =  W(r  =  0),  and 


Ve)  -  pN 


p±sin26  2 


k±4P4tan20 


V(6)  (1  +  k ,  2 p 2 s ^n2 e )  3 


(33) 


Thus,  n(2(Ue>0)  'v  k±4p4  in  the  limit  in  which  k±p  <  1,  and 


n(2u  ,0)  *v  (k,p)  in  the  opposite  case. 


The  emissivity  expressed  in  (32)  and  (33)  admits  relatively 

simple  numerical  analysis,  and  we  display  the  angular  spectra  of  the 

ordinary  and  extraordinary  modes  in  Figs.  1  and  2  by  plotting  <3^(0)  for 

0^6^  tt/2  and  several  choices  of  to  /Q  and  to  p/c.  As  shown  in  the 

e  e  e 

figures,  there  is  a  quadrupole  radiation  pattern  for  both  the  ordinary 
and  extraordinary  modes,  which  is  ..ighly  sensitive  to  the  transverse 
scale  size  of  the  soliton.  Specifically,  for  togp/c  >  1  the  radiation 


is  strongly  beamed  in  the  directions  parallel  and  antiparallel  to  soli- 

ton  propagation,  and  for  «ep/c  ^  1  peak  emission  occurs  for  0  ^  tt/4. 

We  note  that  for  co  p/c  <  1,  the  emissivity  scales  as  n(2w  ,6)  'v  sin220 
e  e 

in  the  limit  of  large  im  /Q  . 2,4  In  Fig.  3,  we  plot  n,(2w  ,0)/n  (2 w  ,0) 

6  G  ">6  ““6 

for  parameters  consistent  with  those  used  in  the  computations  of  Fig.  1. 

We  remark  that  this  quantity  appears  to  be  relatively  insensitive  to  p 

over  the  range  studied  (i.e.,  0.1  <  wep/c  <  10  ) ,  and  we  display  the 

result  for  w  p/c  =  1.  The  principal  results  are  that  (1)  for  m  /fi  <1 
e  e  e 

the  ordinary  mode  tends  to  dominate  the  emission,  but  that  this 
situation  is  reversed  when  ^e/^e  >  1,  and  (2)  that  the  characteristic 
dominance  of  either  mode  is  greatly  enhanced  for  0  ^  40°. 

Finally,  we  observe  that  condition  (31)  is  equivalent  to  the  re- 

2 

quirement  that  wo/nQTe  »  12(ve/c)  ,  and  that  it  is  clear  from  (32)  that 
in  this  limit  n(2<i>e,0)  'v-  Wq. 

In  order  to  treat  the  case  in  which  k+  <  kq,  we  must  rely  on 
wholly  numerical  methods.  To  this  end  we  first  rewrite 


and 


N^sin20(p  ±  sin2§)  „ 

^(0)|  fdx  xJQ(N+xsin0) 


x  each 


T  Trk+cos0 

L  2k 


exp(xc/uiep) 


y 


(36) 


(37) 


The  dependence  of  the  emissivity  at  u  -  we,  and  on  both  0  and  Wq 
(1  2) 

are  contained  in  P+  ’  ,  and  it  is  these  quantities  that  we  evaluate 

here.  Note,  again,  that  the  condition  required  for  the  validity  of  (34) 
is  that  6N+(vg/c)  >  k^  'v  (wQ/noTe) 1/^2 ,  and  that  N2  %  l^n^/nj  'v  WQ/noTe 
o 

while  N  1.  Thus,  it  is  difficult  to  satisfy  this  requirement  for 

ordinarymode  waves,  and  we  restrict  the  analysis  to  consideration  of  the 

extraordinary  mode  emissivity  at  u>  -  (no  restriction  is  necessary  for 

emission  at  2wg) .  It  is  the  existence  of  an  electromagnetic  mode  with 

frequency  u  -  and  index  of  refraction  N  'v  1  which  constitutes  a  major 

distinction  between  magnetized  and  unmagnetized  plasmas. 

We  consider  the  case  of  emission  at  u>  =  cj  first,  and  plot  the 

results  of  the  numerical  integration  of  versus  w0/n0Te  in  Fig.  4. 

It  should  be  noted,  again,  that  the  constraints  on  the  analysis  in  this 

frequency  regime  that  N_(v  /c)  >  k  X  >  X  /p.  In  the  results  presented,  we 

”  6  0  6  6 
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chose  T  =0.1  keV  and  ai  p/c  =  10,  which  imply  that  p  =  715  A  and 
e  e  e 

-3  -2 

1.4x10  ^icA  h  1,4  x  10  .  For  simplicity,  we  have  assumed  that 

o  e 

Tj  *  0  in  the  analysis.  It  is  clear  from  the  figure  that  the  angular 

spectrum  depends  critically  on  both  and  the  soliton  amplitude, 

and  that  no  simple  scaling  law  can  be  found  between  n_(u)e,9)  and 

W  T  .  It  should  be  observed,  however,  that  while  increases  in  plasma 
n  o  e  r 

density  (i.e.,  in  cu^ / )  leave  the  scaling  at  low  levels  of  soliton 
amplitude  relatively  unchanged,  the  scaling  of  the  emissivity  with 
Wo/noTe  and  the  angular  spectrum  of  the  emission  are  substantially 
altered  at  higher  levels  of  Wq. 

In  Figs.  5  and  6  we  plot  the  results  of  a  numerical  integration  of 
(2)  (2) 

P.  and  P  versus  W  /n  T  respecti  .’e !/ .  It  is  clear  from  both 
+  -  o  o  e 

figures  that  (1)  as  Wq  increases  we  recover  the  result  in  (32)  in 
which  n±(2a>e,6)  ~  Wq,  (2)  at  lower  levels  of  Wq  the  emissivity  increases 
faster  than  Wq,  and  (3)  the  angular  spectrum  of  the  emission  is  sensi¬ 
tive  to  the  soliton  amplitude.  We  also  observe  that  while  the 

(2) 

ordinary  mode  emissivity  (i.e.,  P+  )  is  relatively  insensitive  to  the 

plasma  density,  the  extraordinary  mode  emissivity  is  greatly  modified 

in  going  from  w  /SI  =  0.1  to  w  /U  =  10.  This  can  be  explained  by 
e  e  e  e 

noting  that  the  emission  is  in  the  slow  extraordinary  mode  for 

u)  /Q  =0.1  (i.e.,  the  emission  frequency  is  below  the  upper  hybrid 
e  e 

frequency),  and  the  fast  extraordinary  mode  for  =  10.  This  will 

have  severe  consequences  on  the  radiation  observed  from  outside  the 
plasma,  since  slow  extraordinary  mode  waves  cannot  readily  escape  from 
the  plasma  without  tunneling  through  the  upper  hybrid  layer  or  mode 
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coupling  to  the  ordinary  or  fast  extraordinary  modes.  Finally,  we 

remark  that,  as  shown  in  (32),  the  angular  spectrum  of  the  radiation 

and  the  scaling  of  ni(2a)e,6)  with  WQ  should  also  be  sensitive  to  the 

transverse  dimension  of  the  soliton;  however,  it  is  beyond  the  scope 

of  this  work  to  treat  this  scaling  in  the  regime  in  which  k .  >  <  . 

x  o 

IV.  SUMMARY  AND  DISCUSSION 

In  this  work,  expressions  have  been  derived  for  the  radiation  of 

an  arbitrary  three  dimensional  Langmuir  wave  packet  at  u  =  in  and  2u> 

e  e 

in  a  uniformly  magnetized  plasma.  The  analysis  of  the  radiation  at  the 

plasma  frequency  has  been  limited  to  the  regime  in  which  the  radiation 

wavelength  is  much  less  than  the  scale  length  of  the  soliton  in  the 

interest  of  deriving  an  analytic  expression  for  the  emissivity,  which 

2  2 

imposes  the  requirement  that  12(ve/c)  N±  WQ/noTe.  Since  the  only 

electromagnetic  mode  in  a  field-free  plasma  which  frequency  w  =  ioe  has 

2 

an  index  of  refraction  N  ~  W  /n  T  ,  this  condition  imposes  a  severe 

o  o  e 

restriction  on  the  present  analysis  to  that  of  a  very  hot  plasma.  How¬ 
ever,  the  presence  of  an  ambient  magnetic  field  introduces  an  additional 
mode  with  a  mixed  electrostatic/electromagnetic  polarization  (i.e.,  the 
extraordinary  mode)  having  an  index  of  refraction  of  the  order  of  unity 
in  the  vicinity  of  the  plasma  frequency,  and  which  presents  no 
such  severe  restriction.  While  these  waves  cannot  readily  escape  from 
the  plasma  (unless  some  means  of  tunnelling  through  the  upper  hybrid 
layer  or  mode  conversion  to  the  ordinary  or  fast  extraordinary  modes 
is  possible)  and  should  not  be  an  important  characteristic  of  radiation 
from  astrophysical  plasmas,  study  of  this  radiation  mode  may  be 
important  in  laboratory  plasmas.5-9 


In  order  to  determine  relatively  simple  expressions  for  the 
radiation  emissivity  and,  thereby,  to  determine  the  angular  spectrum  of 
the  emission  as  well  as  the  scaling  of  the  radiated  power  with  soliton 
amplitude,  the  specific  case  of  one  dimensional  Langmuir  solitons  has 
been  studied  in  some  depth.  In  this  limiting  regime,  it  is  shown  that 
the  angular  spectrum  is  sensitive  to  both  the  soliton  amplitude  and  to 
the  transverse  scale  size  of  the  soliton.  While  no  simple  scaling  law 
between  the  emissivity  and  the  soliton  amplitude  is  readily  apparent 
for  w  =  <i>e,  it  is  clear  that  for  oo  =  2we  the  emissivity  is  linearly 
proportional  to  the  soliton  amplitude  when  W^  exceeds  a  certain  thresh¬ 
old  which  depends  on  the  plasma  density,  the  ambient  magnetic  field, 
and  the  angle  of  propagation  of  the  radiation.  The  -fmmediate  signifi¬ 
cance  of  this  result  is  to  the  scaling  of  the  second  harmonic 
radiation  in  type  III  solar  bursts.2 
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APPENDIX:  THE  DISPERSION  TENSOR  AND  SOURCE  CURRENT 


The  fluctuation  fields  giving  rise  to  the  emission  are  of  high 

order  in  VA,  and  in  order  to  treat  emission  at  w  and  2ui  we  must  solve 

e  e 

the  equations 


If 


—  <5E  -  v 

m  <b  'V 

e 


(1).  yv(l) 


Q  e 
e  'uz 


6v 

% 


(Al) 


x  6B  «  - 
-b  c 


3_ 

3t 


6E 

% 


[<*.• 


+6n)6v  + 

'b 


(n 


(1) 


6n)v 

a. 


(A2) 


V  x  SE 


-  i  6B 

c  3t  ^ 


(A3) 


where  6E  and  6B  are  the  radiation  fields,  6v  is  the  high  order  velocity 

<b  % 

fluctuation,  and  n^  and  v^  are  the  first  order  density  and  velocity 

'Xj 

fluctuations.  Note  that  Sn,  which  describes  the  caviton  structure,  is 
itself  of  second  order  in  V<t> .  Thus,  the  term  in  6n  Sv  is  of  at  least 
fourth  order  in  Vcj>  and  gives  rise  to  the  turbulent  shift  in  the  plasma 
frequency.  In  addition,  the  term  in  dnv^1^  is  of  third  order  in  V<j> 
and  is  responsible  for  the  oscillatory  current  at  we.  In  contrast, 
the  terms  giving  rise  to  emission  at  2oig  are  of  second  order  in  V<J>. 

Eliminating  6B  from  this  system  of  equations,  we  find  after  some 
straightforward  manipulations  that 


(k  k  -  k2I)  +  I  - 

%  %  fa  fa 


•6E(k,ui) 


m 


e 


e 


o 


<zM 


,01 


Airie 


(dn 


<Sv). 

'o  k,ui 
% 


_  ITe  / 

oi2  \3t 


(n 


(1) 


<5n)v 


0)\ 
h •• 


(A4) 
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COSO)  t 
e 


(1)  _  _JL 


m  0) 

e  e 


v) 


+  ^  2-f  "2  <|2  x  V*)sinu»et 

e  e 


(All) 


Substitution  of  (A10)  and  (All)  into  (A7)  reproduces  Eqs.  (10)  and  (11), 
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W0/n0T, 


W0/n0T. 

Fig.  6  —  Plots  of  Pl2>  versus  W0/n0Te  for  Te  =  .1  keV,  coep/c  =  10,  and 
(a)  we/ne  *  0.1,  and  (b)  we/Sle  =  10 


